SUMMARY 


ii} from the two-phonon processes to the intensity of diffuse X-ray scattering 
has been accurately calculated in the case of aluminum for three different directions 
scattering vector x=k—k,, where k, is the wave vector of the incident photon and k that 


vibrations in the harmonic approximation. The only further approximation introduced 

is the assumption that the temperature is so high that the mean energy of every lattice vibra- 
vid tion may be put equal to kT. Comparisons are made between this accurate calculation 
_ and different approximate calculations. This investigation is of interest for the determination 


fuse X-ray scattering. ' < 


bes a 
_ Measurements of the intensities of diffuse X-ray scattering by a single crystal 
make it possible to determine the frequency spectrum of the crystal. For the general 


the vibrating atoms only in the harmonic approximation, the intensity of the 
_ thermal diffuse scattering can be expressed as a sum of terms with the following 
_ interpretation. The first term represents the contribution from scattering processes 
- involving only one elastic wave (one-phonon processes), the second term involves 
two such waves (two-phonon processes), and so on. In general only the first two 
or three terms contribute appreciably to the intensity. The contribution from the 
_ one-phonon processes is simply related to the frequency spectrum of the crystal. 
_ Experimentally this intensity contribution is obtained from the value of the total 
- measured intensity by subtracting the Compton-scattering and the two- and three- 
_ phonon scattering, effects from higher-phonon processes being negligible [2]. The 

terms which are to be subtracted are sometimes comparable to or even larger than 

the contribution from the one-phonon scattering processes [2]. 

The two-phonon scattering has until now been calculated approximately in various 
ways. Olmer [3, 4], in the case of aluminum, used the simple Debye model, i.e., 
assumed the crystal to be an ideal isotropic continuum, all the elastic waves to 
have equal mean velocities, and the dispersion curve! to be linear. Walker [2], 
also dealing with aluminum, assumed the crystal to be elastically isotropic but 
distinguished between longitudinal and transverse waves and introduced a sine- 


1 The dispersion curve is obtained by plotting the frequency as a function of the wave 
number. 


per atom due to the two-phonon scattering ist a 
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where (33) is Thomson’s expression for the differential scattering cross 
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section of a free electron, and I(x) is the integral. 
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Here the vectors q and q’ are the wave vectors of the two phonons participating 
in the scattering process. The vector x is equal to the difference between the 
wave vector k of the scattered photon and the wave vector k, of the incident 
photon. The vectors q and q’ are to be chosen to lie within the first Brillouin 
zone and to fulfill the condition 
q+q' =x-Tt (3) 

where t is a reciprocal lattice vector.2 For each t the integration in (2) is there- 
fore to be performed over the volume V, common to the first Brillouin zones 
centered on the scattering point (i.e., the endpoint of x) and on the reciprocal 
lattice point (i.e., the endpoint of +t). When the sum over all possible t is per- 
formed, these volumes exactly fill the first Brillouin zone centered on the scattering 
point. The other notations appearing in (1) and (2) are defined as follows: 

/, is the atomic scattering factor 

e *™ is the Debye-Waller temperature factor 

m is the atomic mass 

V is the volume of the direct unit cell 

k is Boltzmann’s constant 

T is the absolute temperature 

e, 1s the unit vector in the direction of x 

«;(q), where 7=1, 2, 3, are the frequencies of the three elastic waves correspond- 
ing to the phonon wave vector q 

C;(q), where j=1, 2, 3, are the amplitude vectors of the three elastic waves cor- 
responding to the phonon wave vector q. These amplitude vectors are normalized 

* For the general formulae we refer to Waller [1]. We remark that Walker’s [2] formula (4) is 
the same as our formula (1). 

* Our definition of the reciprocal lattice differs from the usual one. Cf. the text to isa ie 
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ig. 1. This figure illustrates the relations in the reciprocal space of the aluminum lattice. The 
reciprocal lattice, which is body-centered cubic, has lattice points at the corners and at the 
of the cube, the edge-length of which is 4/a. (According to the usual definition of the 
ocal lattice the edge-length is 2/a, but we have introduced a factor 2 for convenience.) 
polyeder represents the first Brillouin zone. 


such that 
Res 2 Che Coy = Oxy (4) 


where 6,, is the Kronecker symbol. 
Starting from the dynamical equations for the vibrations of a monatomic lattice 


2 (Dy — Oxy 3) Cy =0 (5) 


where D,, are quantities which depend on the interatomic forces, we transform 
(2) into a form which is more convenient for our purpose [1]. Multiplying (5) by 
Bis 
ie 
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summing over 7, and using (4), we get 
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one gets ie: 
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where A (q) denotes the determinant having the elements D,, (q), and A,, (q) denotes — 
the cofactor belonging to the element D,, (q). Substituting (7) into (2), we obtain 


| 3o=to daw | [Reet Aaa) | 


ices AWA G=«) s e 


First Brillouin 
zone 


where use has been made of the facts that A(q) and A,,(q) are periodic in the 
reciprocal lattice. This fact is immediately clear, for the elements D,, (q) are known 
to have this periodicity. 

We shall now consider the formula (8) for the case of aluminum, which has a 
face-centered cubic lattice. The reciprocal lattice is body-centered cubic, and the 


volume of integration in (8), which is the first Brillouin zone shown in Fig. 1, 
is limited by the planes 


Wr 4 3 
dy re +t951+ \o,)+1¢. = 
qz 


where q,, q,, and q, are the cartesian components of the wave vector q, and a is 
the lattice constant. Considering only the forces between nearest and next nearest 
neighbours and assuming the forces of the latter kind to be central forces, one has 
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Fig. 2 c. 


‘Fig. 2. This figure shows plots of I (%), expressed in sect cm~*, as functions of | * |, expressed in 

_ units of 2/a, where a is the lattice constant. The figures 2a, 2b, and 2c refer to the cases that 
the scattering vector % lies along the directions (100), (110), (111), respectively. Along these three 
directions I (%) is symmetric with respect to the middlepoint of the line joining two neighbouring 
lattice points and is also periodic in the reciprocal lattice. The curves A, B,C, and D, which are 
drawn up to the symmetry point, have been obtained as follows: 

A. Numerical evaluation of the integral (8). 

B. Use of Walker’s method for evaluating the integral (2). _ 

C. Use of the simple Debye model for evaluating the integral (2). 

D. Use of a simplification of Walker’s method for evaluating the integral (2). 
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04, 0%» By, Y being interatomic force constants, the values of which are taken 
from Walker [2] (cf. his table I, model -B). ae 


In Fig. 2 the integral J (x) is plotted against |x| for the case of aluminum and 
for three different directions of x. Fig. 2a refers to the direction (100), Fig. 2b 
to the direction (110), and Fig. 2¢ to the direction (111). For each direction the 
integral I(x) has been calculated in four different approximations, and the results 
are represented by the curves A, B, C, and D. These curves have been obtained as 
follows: 

The curves A. These curves are obtained by computing the integral (8) numerically 
using an electronic computer. 

The curves B. These curves are obtained according to a method which has been 
introduced by Walker [2] and in which the following approximations are made: 

1°. The crystal is assumed to be an isotropic continuum, and consequently the 
phonon waves are considered to be purely longitudinal or purely transverse. This 
property of the waves is true for small values of |q| but certainly not for large 
values. 

2°. The dispersion curves, i.e., the curves giving the frequency as function of 
the wave vector, are assumed to be sine-shaped, i.e., 


2Om ef ae 
oy(a)=v, Hsin (Zl), 0<|a|<an (9) 
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v; being the sound velocity of the longitudinal (j=1) and transverse (j= 2,3) 
waves, respectively, and q,, being the maximum value of | q| as defined in 3° below. 
From the actual dispersion curve it is seen that this approximation is not good 
in all directions (cf. the figures 2, 3, and 4 in Walker’s paper [2]). 

3°. The actual volume of integration, i.e., the first Brillouin zone (Fig. 1) is 
replaced by a sphere of the same volume, the radius of which is denoted by qn. 

Walker’s resulting formulae are given in the appendix in his paper [2] (cf. his 
eqs. (10) and his table II). Our curves B are calculated from the right-hand 
side of Walker’s eq. (10), which we have, however, multiplied by the factor 


3 1\2 
kied (> ~3} in order to get values corresponding to our J (x). Comparing the curves 
9 Om j=l Vj 
A and B, we see that Walker’s approximate method gives astonishingly good results. 
The curves €. These curves are obtained by using the simple Debye model for 
evaluating the integral I(x). This means that the following approximations are 
made: 
i The crystal is assumed to be elastically isotropic, but no distinction is made 


between longitudinal and transverse waves, i.e., [e,-C;(q)? is put equal to 1/3 
for all. j= 1, 2, 3. 
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2°. The relation between frequency and wave vector is approximated by 


o;() =I 4] (10) 
3°. The actual volume of integration is replaced by a sphere of the saine volume. 


If these approximations are introduced, one can perform two of the three in- 
tegrations in (2) analytically, getting 


a / 3 1\2 gnl|x—t| . 
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From this formula the curves C are obtained. It is seen that they give much too 
low values of I(x), especially for large values of |x|. 
The curves D. These curves are obtained by modifying Walker’s method (from which 
the curves B were obtained) such that we replace Walker’s sine-shaped dispersion 
curves (9) by the linear dispersion curves (10) appearing in the simple Debye model. 
The integral (x) can then be evaluated analytically, except for one term which 
contains the integral (11) appearing in the simple Debye model. The curves D give 
still lower values of J (x) than the curves C. The simple Debye model (curves C) is 
thus seen to give better values of J (x) than the (in a sense more detailed) model 
used to obtain the curves D. The use of the linear approximation (10) for the 
dispersion curves causes w;(q) to be larger than its actual value, especially for large 

values of |q|. In the linear approximation (10) the denominator [a; (q) «;(q')/ of 
the integrand on the right-hand side of (2) is therefore too large, and it is clear 
that the curve D must lie below the corresponding curve B. In the approximation 

corresponding to the curves C the too large values of the denominator of the 
integrand are somewhat compensated by too large average values of the numerator 
[C; (q) - ex]? [C; (q’)- e.]?, which appear when each factor of this numerator is replaced 
by 1/3. When the actual angular dependence of the numerator is taken into account 
in detail, as is done in the calculation of the curves D, no such compensation 
is introduced. This is the reason why the curves C are better than the curves D. 
Although the approximation corresponding to the curves D gives too low values, 
these curves have approximately the right shape. It is in fact possible to get curves 
which are similar to the curves A by multiplying the values of J (x), according 
to the curves D, by a constant factor of about 2.3. This corresponds to using, 
instead of (10), the formula 


(OF (q) ~ 0.8 2;| q| (12) 
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